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We present simulation results for lattice QCD with chiral fermions in small volumes, 
■ where the e-expansion of chiral perturbation theory (%PT) applies. Our data for the 
low lying Dirac eigenvalues, as well as mesonic correlation functions, are in agreement 
with analytical predictions. This allows us to extract values for the leading Low Energy 
Constants F n and E. 

1. The e-regime of QCD 

For QCD at low energy and zero quark masses one assumes the chiral symmetry to 
break spontaneously, SU(N f ) L ® SU(N f ) R -> SU(N f ) L+R . This yields Nf - 1 Nambu- 
Goldstone bosons, which pick up a small mass when we proceed to light quark masses m q 
(we assume the same m q for all flavours involved). These bosons are identified with the 
c-j \ light mesons, in particular the pions (and often also the kaons and if). 

%PT is a powerful description for them by means of an effective chiral Lagrangian, 

C eS [U] = -f Tr[d,Ud^] - ^ S Tr[U + rf\ + . . . , 

which deals with fields U in the coset space SU(Nf) of the chiral symmetry breaking. 
One writes down the terms which fulfil the symmetry requirements and classifies them 
according to low energy counting rules. The coefficients to these terms are the Low Energy 
Constants (LEC), hence F n and S are the leading LEG 

One usually considers this effective theory in an infinite or very large volume, which 
allows for an expansion in the meson momenta (p-expansion) . However, one can also 
consider the opposite situation — a small volume with linear size L <C 1/ — which is 
called the e-regime. Here the p-expansion fails due to the importance of the zero modes. 
Fortunately the path integral for the latter can be performed with collective variables [ 
Then an expansion in the low (non-zero) momenta and the meson mass is possible 
(e-expansion). 

It is in fact motivated to study the e-regime, although it is not a physical situation to 
squeeze pions into such a tiny box. The point is that even in a small volume values for the 
LEC in infinite volume can be determined — exactly by analysing the finite size effects. 

The LEC enter the effective Lagrangian as free parameters, hence they can only be 
determined from QCD as the fundamental theory (or from experiment). This is a goal 
for lattice simulations, which are clearly more convenient in small volumes. 
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However, simulations in the e-regime face technical difficulties, which could only be 
overcome in the recent years. One has to keep track of the chiral symmetry and deal with 
very light pions. Moreover, we need a clean definition for the topology, since expectation 
values should often be measured in distinct topological sectors These properties 

are provided by Ginsparg- Wilson fermions, which have a lattice modified, exact chiral 
symmetry and exact zero modes, so that the topological charge can be defined as the 
index v (i.e. the difference between the number of zero modes with positive and negative 
chirality) [Hj. The simplest realization of a Ginsparg- Wilson fermion uses the overlap 
Dirac operator [E]. 

Although algorithmic tools for the e-regime have been worked out [E], such simulations 
are still computationally expensive: for the time being one has to use the quenched 
approximation — which leads to logarithmic finite size effects on the LEC [ [7] . 

2. The Dirac spectrum 

Chiral Random Matrix Theory (RMT) simplifies QCD to a Gaussian distribution of 
fermion matrix elements, and still manages to capture some aspects of QCD. One of 
its numerous assumptions is that the energy should be well below the Thouless energy 

-^Thouless = F%/(\/VY). 

In particular, Ref. [IHj presents explicit predictions for the microscopic spectral density 

oo 

= p} (z), z = XEV, A : eigenvalue of the Dirac operator. 

The density in each topological sector \v\ is further decomposed into densities of the 
leading individual eigenvalues, px"(z) = Xm=i,2,3... Pn{z). Using the overlap operator 
at m q = we have measured P\\z) at \v\ = 0, 1 and 2 in different volumes [E]. In 
Figure n we compare our data for the cumulative density with the RMT predictions. The 
comparison fails in the small volume V ~ (0.98 fm) 4 , but it works well in V ~ (1.23 fm) 4 . 
In this fit, S appears as the only free parameter. The successful fit in V ~ (1.23 fm) 4 
implies a value of £ ~ (253 MeV) 3 , which agrees well with the literature. 

The microscopic RMT predictions work in such volumes up to the second eigenvalue, 
then the spectral density starts to grow cubically (bulk behaviour). This transition can 
be viewed as an effective bound for the chiral RMT, thus taking the role of Erhouiess- Our 
observations were confirmed by subsequent studies, which included larger volumes [HUj. 

3. Axial vector correlation functions 

Since QCD simulations with chiral fermions are restricted to the quenched approxi- 
mation at present, also quenched %PT was worked out [ as a basis of predictions 
that the lattice data can be confronted with. Regarding the mesonic correlation func- 
tions in leading order, the vector correlator vanishes (actually to all orders), and the 
scalar and pseudoscalar correlators have the drawback that they involve additional LEC, 
which are specific to quenching. Hence we considered the correlator of the axial current 
A pit) = ^g^(t,x)j5j^(t,x), where the leading order only involves F n and S. The e- 
regime prediction for {A (t)Ao(0))\ u \ in a volume L 3 x T is a parabola in t with a minimum 
at T/2, and with an additive constant F^/T (this is in contrast to the usual cosh be- 
haviour in large volumes). The curvature in the minimum of the parabola is oc S. We now 
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Figure 1. The cumulative density of the first non-zero eigenvalue in the sectors \u\ = 0, 1 
and 2. The comparison of our lattice data to the chiral RMT predictions (curves) fails in 
V = (0.98 fm) 4 (on the left), but it is successful in V = (1.23 fm) 4 (on the right). In the 
latter case, the fit yields £ = (253 MeV) 3 . 



simulated overlap fermions at m q = 21.3 MeV 



A first experiment with L ~ 0.98 fm 



failed to reproduce the prediction. This is exactly consistent with the observation for the 
Dirac eigenvalues described before: in both cases we badly need L > 1.1 fm. 

We also observed that the Monte Carlo history at v = is plagued by spikes, so that a 
huge statistics is needed. These spikes occur at the configurations which have a very small 
(non-zero) Dirac eigenvalue. As Figure ^ shows, this probability is strongly suppressed 
for \v\ > 0, and indeed these histories are smoother. 



So we measured in a volume V = (1.12 fm) 4 at 



1. Figure|21 (on the left) shows 



that we found a nice agreement with the predicted parabola, which is plotted for E = 
and for £ = (250 MeV) 3 . Since these curves can hardly be distinguished, we see that the 
sensitivity of the curvature is by far too small for a determination of £ in this way. On 
the other hand, this result provides a stable value of F n = (86.7 ±4.0) MeV. The plot on 
the right-hand-side of Figure El shows the value we obtained as a function of the fitting 
range around the minimum. The result that we obtained for \v\ — 2 (with lower statistics) 
agrees within the errors. The renormalisation due to the factor drives our result for 
F n up to a value around 130 MeV, but it is well-known that quenched results for the pion 
decay constant tend to be significantly above the experimental value of ~ 93 MeV resp. 
w 86 MeV in the chiral limit [IT3]. 
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Figure 2. The axial correlation function (on the left) and the values of F w obtained from 
fits over tf points around the centre t = 6 (on the right). 



4 



4. Conclusions 

We have performed two pilot studies about the feasibility to determine LEC of 
from simulations with chiral lattice fermions in the e-regime. Our results show that such 
lattice data can be fitted to the prediction by chiral RMT and quenched However, we 
also found crucial conditions for this to work: the spatial box length must exceed a lower 
bound, L>1.1 fm. Moreover, correlators should not be measured in the topologically 
trivial sector. 

can be extracted from the axial vector correlator, and E from the densities of the 
lowest eigenvalues of the Dirac operator. Thus our pilot studies confirm that simulations 
in the e-regime do have the potential to evaluate LEC from first principles. One could 
also deal solely with the zero-mode contributions to the mesonic correlation functions [ 
d. 

Currently a non-standard lattice gauge action is being tested, which preserves the topo- 
logical charge over long intervals of the Monte Carlo history (far beyond the autocorrela- 
tion time of other observables) [ . This would allow us to collect many configurations 
for measurements in a specific topological sector, which is tedious with the standard gauge 
action. On the fermionic side, we are now including an improved overlap operator [El] in 
our investigation of the e-regime. At last we add that we also have results in the p-regime 
[I17j. for instance for the direct measurement of F n . 

This work was supported by the Deutsche Forschungsgemeinschaft through SFB/TR9- 
03. Most computations were performed at the ZIB/HLRN in Berlin. 
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